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Knots in a Spinor Bose-Einstein Condensate 
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We show that knots of spin textures can be created in the polar phase of a spin-1 Bose-Einstein 
condensate, and discuss experimental schemes for their generation and probe, together with their 
lifetime. 
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A decade ago, Faddeev and Niemi suggested that knots 
might exist as stable solitons in a three-dimensional clas- 
sical field theory, thus opening up a way to investigate 
physical properties of knot-like structures |l[ . They fur- 
ther proposed Q that their model can be interpreted 
as the low-energy limit of the Yang-Mills theory, where 
knots are suggested as a natural candidate for describing 
glueballs - massive particles made of gluons. 

In cosmology, topological defects are considered to be 
important for understanding the large-scale structure of 
our universe 0]. Although recent measurements of the 
cosmic microwave background (CMB) have shown that 
topological defects are not the dominant source of CMB 
anisotropics, the search for topological objects in the uni- 
verse, such as cosmic strings, continues to be actively 
conducted. Recently, it has been suggested that a cosmic 
texture, classified by the homotopy group of tt 3 (S 3 ) = Z 
Q, generates cold and hot spots of CMB Q. This tex- 
ture is a spherical or point-like object that is unstable 
against shrinkage according to the scaling argument. In 
cosmology, the instability is favored for evading cosmo- 
logical problems such as a monopole problem. Three- 
dimensional skyrmions [B| and Shankar monopoles [E @] 
are topological objects belonging to the same homotopy 
group. On the other hand, knots, which belong to a dis- 
tinct homotopy group, tts(S ) = Z, have thus far been 
ignored by cosmologists; however, they would be a po- 
tential candidate for topological solitons in our universe. 

Knots are unique topological objects characterized by 
a linking number or a Hopf invariant as discussed in 
a seminal paper on superfluid 3 He Q. Other famil- 
iar topological objects, such as vortices JE |E EE Ell ], 
monopoles [E EH, and skyrmions [E @, U are charac- 
terized by winding numbers, which have recently been 
discussed in relation to spinor Bose-Einstein condensates 
(BECs) [H El EE EE EE EE El- However, little is 
known about how to create such knots experimentally. 
In this Letter, we point out that spinor BECs offer an 
ideal testing ground for investigating the dynamic cre- 
ation and destruction of knots. We also show that knots 
can be imprinted in an atomic BEC using conventional 
magnetic field configurations. 



We consider a BEC of spin-1 atoms with mass M that 
are trapped in an optical potential V(x). The energy 
functional for a BEC at zero magnetic field is given by 



E= d 



h 2 

— y iv^ n 

2M ^ 1 v 



t/ , go 2 

VP+-P 



9i 21*1: 
yP l/l 



(1) 



where ip m (x) is an order parameter of the BEC 
in a magnetic sublevel m = 1,0 or —1 at posi- 
tion x, and p{x) = E m lV>m(aO| 2 and p(x)f(x) = 
Emm'C^l^im'^m'^) are the number density and 
spin density, respectively. Here F = (F x ,F y ,F z ) is 
a vector of spin-1 matrices. The strength of the in- 
teraction is given by go = Airh 2 (2a2 + ao)/(3M) and 
gi = 4irh 2 (a2 — ao)/(3M), where 05 is the s-wave scat- 
tering length for two colliding atoms with total spin S. 
The ground state is polar for g\ > and ferromagnetic 
for 5 i < [IE [21|. 

The order parameter for the polar phase can be de- 
scribed by the superfluid phase $ and unit vector field 
n in spin space, whose components are given in terms of 
1pm as 
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This order parameter is invariant under an arbitrary ro- 
tation about n = (n x ,n y ,n z ), i.e., exp[— in ■ F^]^ = 
where 7 is an arbitrary real number. It is also invari- 
ant under simultaneous transformations $ — > iJ + it and 
n — > — n. The order parameter manifold for the polar 
phase is therefore given by M = (U(l) x 5 2 )/Z 2 
where U(l) denotes the manifold of the superfluid phase 
i?, and S 2 is a two-dimensional sphere whose point spec- 
ifies the direction of n. 

Knots are characterized by mappings from a three- 
dimensional sphere S* 3 to S 2 . The S 3 domain is pre- 
pared by imposing a boundary condition that ^ takes 
on the same value in every direction at spatial infin- 
ity, so that the medium is compactified into S 3 . Since 
neither U(l) nor Z2 symmetry contributes to homotopy 



groups in spaces higher than one dimension, we have 
7r 3 (M) = ^(S 2 ) = Z. The associated integer topological 
charge Q is known as the Hopf charge which is given by 
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where Ti. 
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n ■ (diti x djfi) [lj. Note that 



the domain (a;) is three-dimensional, while the target 
space (n) is two-dimensional. Consequently, the preim- 
age of a point on the target S 2 constitutes a closed loop 
in S 3 , and the Hopf charge is interpreted as the linking 
number of these loops: If the n field has Hopf charge 
Q, two loops corresponding to the preimages of any two 
distinct points on the target S 2 will be linked Q times. 

Figure 1 (a) illustrates the n field of a knot created in 
a uniform BEC with the linking number of 1. Here, the n 
field is expressed as n(x) = exp —ia(r)v(6, </>) • F 

subject to the boundary condition = (0, 0, 1) T at spa- 
tial infinity, where x = (r sin 9 cos <f>, r sin 9 sin 0, r cos 9) 
and F is a vector of spin-1 matrices in the Cartesian rep- 
resentation, i.e., Fj - 



mapping v(p, q. 
and we choose 



F j = ~ ie jki {j,k,l = x,y,z). A given 
i : S 2 — > S 2 determines the charge Q, 
v = (— sm9 cos 0, — sin 6* sin (f>, cos 6*) in 
Fig. 1 (a). The radial profile function a(r) is a monoton- 
ically decreasing function of r, subject to the boundary 
conditions a(0) = 2ir and a(oo) = 0. Here, we consider 
a(r) = 2tv[1 — tanh(r/£knot)]> where £knot is a character- 
istic size of the knot. Although there is no singularity in 
the texture, it is impossible to wind it off to a uniform 
configuration because this texture has a nonzero Hopf 
charge of 1 . The left part of Fig. 1 (b) describes the or- 
der parameter for the to = — 1 component in real space, 
where we plot the isopycnic surface of the density and 
the color on the surface represents the phase. Here, n on 
the torus is almost perpendicular to and the phase 
of ip-i is given by arg(i/>_i) = arctan(n !/ /ri x ). The core 
of the knot (n = — rioo) is depicted as a white tube in 
Fig. 1(b). On the right side of the figure, we show the 
extracted preimage for n = (0, 0, — 1) T (white tube) and 
that for n = (1, 0, 0) T (black tube). 

The torus shape of the m = — 1 component appears 
as a double-ring pattern in the cross-sectional plane at 
z — 0, as shown in Fig. 1 (c), where the density distribu- 
tions of the to = ±1 (left) and to = (right) components 
are shown in the gray scale. The distributions of to = ±1 
components overlap completely; therefore the system re- 
mains unmagnetized. This double-ring pattern can serve 
as an experimental signature of a knot and it should be 
probed by performing the Stern-Gerlach experiment on 
the BEC that is sliced at z = 0. 

Next, we show that knots can be created by manipu- 
lating an external magnetic field. In the presence of an 
external magnetic field, the time-dependent phase dif- 
ferences between different spin components are induced 
because of the linear Zeeman effect, which causes the 
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FIG. 1: (Color) (a) Three-dimensional configuration of the 
n field of a knot with Hopf charge Q = 1, expressed as 

n = e- lQ(r),i ' i? (0,0, 1) T , where a(r) = 2tt[1 - tanh(r/&nat)] 
and v = (— sin 9 cos <f>, — sin 9 sin <j>, cos 9). Note that n = 
(n x ,n y ,n z ) T is related to the order parameter St through 
Eq. (J5J. In the figure, only n on the x,y, and z axes is 
shown as arrows. The dashed line traces the point n = 
(0, 0, — 1) T , which is a circle on the xy plane with the ra- 
dius [tanh _1 (l/2)] £knot — 0.55 £knot- The color on the ar- 
rows represents the value of n z (see the scale), (b) The torus 
shape on the left side of the figure shows the isopycnic sur- 
face of |i/'-i| 2 = 0.47, which is equivalent to the region where 
\n z \ < 0.24. The color on the surface shows the phase of tp-i, 
that is equivalent to arctan(n y /n x ). The white and black 
tubes are the preimages of n ~ (0, 0, — 1) T (n z < —0.95) 
and n ~ (1,0, 0) T (n x > 0.95), respectively, which are repro- 
duced on the right side of the figure. The two tubes cross 
once (linking number 1), which is consistent with Q = 1. (c) 
Cross sections of the density of the m = ±1 (left) and m — 
(right) components in the xy plane. The size of each panel is 
4£knot x 4£knot- The torus shape of the m = ±1 components 
appears as double rings in the cross section, which can serve 
as the signature of a knot. 



Larmor precession of n, while n tends to become paral- 
lel to the magnetic field because of the quadratic Zee- 
man effect. Suppose that we prepare a BEC in the 
to = state, i.e., n = (0,0, 1) T , by applying a uni- 
form magnetic field Bq in the z direction. Then, we 
suddenly turn off Bq and switch on a magnetic field, 



B(x) = b(r) (sin 8 cos(4> + cj) ), sin 9 sin((f> + (j) ), cos 9) T , 
where b(r) is an arbitrary function and <po is a real pa- 
rameter. This magnetic field configuration is quadrupo- 
lar if <pQ = 7r and monopolar if <j)Q = 0. In what follows, 
we shall consider the case of <j>o = it and b(r) — b'r, un- 
less otherwise specified. Because of the linear Zeeman 
effect, n starts rotating around the local magnetic field 
asn = exp[ — % (li B /2)B(x)-Ft/h}{0,0,l) T , and therefore 
the n field winds as a function of t. 

Figure 2 shows the dynamics of the creation and de- 
struction of knots in a spherical trap subject to the 
quadrupole field. Figures 2 (a)-(d) show the snapshots 
of (left) and preimages of n = (0,0, —1) T and 

n = (1,0, 0) T (right) at (a) t = 0.5T L , (b) 1.1T L , (c) 
2.2T L and (d) 3.3T L , where T L = 2h/{n B b'R TF ) is the 
period of the Larmor precession at the Thomas-Fermi 
radius -Rtf- The solid circle indicates the periphery of 
the BEC on the xy plane. Knot-like objects enter the 
BEC from its periphery [Fig. 2 (a)], and the number of 
knots increases as the n field winds more and more with 
time [Fig. 2 (b), (c)]. With the passage of further time, 
however, the knot structure is destroyed, as shown in 
Fig. 2 (d). The knots therefore have a finite lifetime, 
which will be discussed later. Figures 2 (e)-(h) show 
cross sections of the density for m = — 1 (top) and m = 
(bottom) components on the xy plane at (e) t — 0.5Tl, 
(f) 1.1T L , (g) 2.2T L , and (h) 3.3T L . We find that as the 
n field winds with time, the number of rings increases. 
This prediction can be tested by the Stern-Gerlach ex- 
periment. 

The knot soliton in the simple nonlinear a model with- 
out a higher derivative term is known to be energetically 
unstable, since the energy of the knot is proportional 
to its size £knot, which can be calculated by integrating 
the kinetic energy ~ Ckifot m the v °l ume of the soliton 
~ Cknot- Knots will therefore shrink and finally disap- 
pear. In the case of an atomic gaseous BEC, however, 
the total energy of the system has to be conserved; there- 
fore, the above-mentioned energetics argument does not 
imply the instability of knots in trapped systems. 

The dominant mechanism for destroying the knot in 
the spinor BEC is the spin current caused by the spatial 
dependence of the n field. The spin current induces local 
magnetization according to 

= -V • jTpi„ = J-! E e ^A^V 2 n A , (4) 

thereby destroying the polar state. The initial polar state 
will spontaneously develop into a biaxial nematic state, 
and eventually result in a fully polarized ferromagnetic 
domain. While n is well-defined in a biaxial nematic 
state as one of the symmetry axes, it is ill-defined in the 
fully-magnetized region. 

Substituting n = e^W^^O, 0, 1) T as used for Fig. 1 
into Eq. (fj|, we analytically calculate the time deriva- 
tive of the local magnetization, whose maximum value is 
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FIG. 2: (Color) Dynamics of the creation and destruction 
of knots in a spherical trap under a quadrupolar magnetic 
field, (a)-(d): The left sides of the figure show the isopycnic 
surface of \ip-i\ 2 /p — 0.4, and the color on the surface shows 
the phase of t/'-i at (a) t = 0.5T L , (b) I.IXl, (b) 2.2T L , and 
(d) 3.3T L , where T L = 2h/(fj, B b' Rtf). The white and black 
tubes on the right are the preimages of n ~ (0, 0, — 1) T (n z < 
—0.95) and n ~ (1,0, 0) T (n x > 0.95), respectively, and the 
solid circles represent the periphery of the BEC [p/p(x = 
0) = 0.01] in the xy plane, (e)-(h) Cross sections of the 
density normalized by p(x = 0) for m = — 1 (top) and m — 
(bottom) components in the xy plane at (e) t = 0.57l, (f) 
1.1T L , (g) 2.2T L , and (h) 3.3T L . 
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FIG. 3: Time evolution of the maximum of the local polar- 
ization in a spherical trap with /?tf = 4.3£knot, starting from 
the knot shown in Fig. 1. The time is measured in units of 



M£ knot /7i. The dashed line shows \f(x) 



£/ r knot (see 



text), which agrees well with the numerical result. The insets 
show |/(a;)| in the xz plane, and the number below each panel 
shows the time elapsed in units of M^ not /h. The size of each 
panel is 1.9£ knot x 3.8£ knot . 



given by \f t f{x,t = 0)| max ~ 14.3ft/(M^ not )- As thc 
magnetization initially increases linearly as a function of 
time (see Fig. 3), we define the lifetime of a knot as 
Tknot = [14-3ft/(M^ not )] _1 . We also numerically calcu- 
late the dynamics of the knot shown in Fig. 1 by solving 
the time-dependent Gross-Pitaevskii equation; the result 
is shown in Fig. 3. In Fig. 3, we plot the maximum 
polarization as a function of time. The dashed line rep- 
resents \f(x)\ max = t/rknot, which agrees well with the 
numerical results. The insets of Fig. 3 show the distribu- 
tion of | /(a?) | on the xz plane. The ferromagnetic domain 
emerges because of the spin current caused by the spatial 
dependence of n, and then expand outward. In the case 
of a 23 Na BEC, the lifetime for a knot with £knot = 10 /im 
is Tknot = 2.5 ms. The lifetime increases with increasing 
the size of knots. 

Finally, we point out that knots in a BEC may 
be used as an experimental signature of a magnetic 
monopole. A magnetic monopole induces the magnetic 
field B = fij (2er 2 )f, which acts in a manner similar to 
the quadrupole field and creates knots. Although b(r) in 
this case diverges at r = 0, it forms a knot-like structure 
on a large scale. For instance, knots expand up to 10 /im 
in the period 43 /is of Larmor precession at r — 10 /im, 
where B = 33 mG. 

In conclusion, we have shown that a spin-1 polar Bosc- 
Einstein condensate can accommodate a knot, which is 
also be shown to be created using a quadrupolar mag- 
netic field. Contrary to knot solitons known in other 
systems, the knots in spinor BECs are immune from en- 
ergetic instability against shrinkage, because the energy 



of the system is conserved; however they are vulnerable 
to destruction caused by spin currents because of the n 
texture. The lifetime of a knot increases in proportion 
to the square of the size of knots and is shown to be 
sufficiently long to be observed in a Stern-Gerlach exper- 
iment. 
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